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C. LHUILLIER, J.-B. FOUET 
We describe four phases found in two-dimensional quantum antiferromagnets. Two 
of them display long range order at T = 0: the Neel state and the Valence Bond 
Crystal. The last two are Spin-Liquids. Properties of these different states are 
shortly described and likely conditions of their occurence outlined. 

1 INTRODUCTION 

Investigating various spin-1/2 systems in two dimensions (2d), by means of exact 
diagonalization (ED) on small samples of spins, we have found, up to now, four 
kinds of antiferromagnetic phases at T=0. Two of them display long range order 
(LRO): Neel or Valence Bond Crystal (VBC) LRO. The others are two different 
Spin-Liquids. The Neel case may be qualified as semi-classical: it is a ground- 
state of the (spin-oo) classical model with an order parameter reduced by quantum 
fluctuations. The others are purely quantum phases. 

The hamiltonians considered are SU{2) invariant, like the Heisenberg hamilto- 
nian and its generalizations as the Ji — J2 model or the Ji — J2 — J3 model and 
the multiple-spin exchange (MSE) model. All these models of magnetism picture 
exchange (or super-exchange) between fermions. The permutation operator of two 
spins reads: 

= 2S, • S, + i (1) 
This yields the Heisenberg hamiltonian: 

H=J2S,-S, (2) 

<v> 

where the sum runs over nearest neighbor pairs of spins on sites i and j , the Ji — J2 
model: 

17 = Ji ^ S, • Sj + J2 Yl • (3) 

<ij> <<ij» 

where the second sum runs over next-nearest neighbor pairs and the Ji — J2 — J3 
model with additional 3th nearest neighbor interrations. We were also motivated 
to study MSE models involving exchange of more than two spins. Among them we 
shall focus on the simplest which brings new physics which involves up to 4-spin 
exchange. Its hamiltonian reads: 

H^J,J2 P'l (P^^>^^ + Pijli) (4) 

<ij> <ijkl> 
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where Pijki stands for the permutation of 4 spins on sites i,j, k, I, P^jli is the inverse 
permutation (note that = ^ij^)- lattices are the square, honeycomb, 

triangular and kagome lattices. 

In order to distinguish between the different phases of such models, from ED 
calculations on small samples (of up to A'^ = 36 spins with present day computers), 
the analysis of the symmetries and scaling properties of the low-lying eigen-levels 
of the spectra is a very efficient approach. Typical spectra of systems in the four 
phases identified are displayed in Fig |l| as a function of the eigenvalues S{S + 1) of 
the square of the total spin (a good quantum number as the hamiltonians are 
SU{2) invariant). A superficial analysis allows an immediate detection of obvious 
distinctive features: 

• Fig |l|a: Neel spectrum, the lowest eigen-levels in each S sector form a separate 
set with an energy increasing as S{S -\- 1) (the quantum origin of the Landau 
free energy per spin / increasing with the magnetization m — S/{N/2) as 
/(m) = /(O) -|- m^/2x + • • • for small m, where x is the spin susceptibility). 

• Fig 0b and |l|c: VBC or type I Spin-Liquid spectra with a gap A between 
different sector of spins {f{m) — /(O) -I- Am/2 -I- • • •) 

• Fig |l]d: type II Spin-Liquid spectrum with a spin-gap and a gapless continuum 
of singlets excitations 

We now briefly sketch the spectral properties of the four different phases and outline 
the groups of systems where they have been found. 



2 Semi-classical Neel phase 



Neel LRO breaks SU{2) and various space symmetries at the thermodynamic limit. 
This is revealed in the spectra by a basic set of low-lying eigen-levels (states on 
dashed line in Fig|l^), called QDJS for Quasi Degenerate Joint States in Refs. 
with energies ~ S{S + l)/2Nx at least for S values up to \/N. They collapse to- 
gether to the absolute ground-state as and separate from the softest magnons 
(dotted line in Fig ^) which decrease as 1/ V^. The QDJS have different S values 
which enable breaking of SU (2) symmetry. The number of QDJS per S value is 
caracteristic: it is one for collinear LRO and equal to the number of way of adding 
p spins of size N/2p for noncoUinear LRO with p sublattices. The QDJS belong 
to different irreducible representations of the space group of the lattice which al- 
low and reveal the space symmetry breaking occuring. For instance, in case of 
the Heisenberg model on the square lattice of Fig pi, where translation symmetry 
is broken, the QDJS have wave vectors k = (0, 0) for even S and k = (tt, tt) for 
odd S values and are invariant under all operations of the point group (taking the 
origine on a site), as expected for this collinear Neel LRO. ED and/or Monte-Carlo 
calculations have shown that the Heisepiherg model on the 2d bipartite lattices iru 
vestigated up to now: square (see Refs.QLl'Q and references therein), honeycomb ,LrEl 
1/5 depleted squareO- • •, displays coUinear Neel LRO. In these unfrustrated sys- 
tems, Monte-Carlo methods can deal now with very large systems which enable 
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Figure 1. Eigen-energies vs S{S + 1) of TV = 36 samples, (a) Spectrum of the Heisenberg model 
on the square lattice, characteristic of coUinear Neel LRO: the lowest energies in each spin sector 
S increase as S{S + 1) with S (states on the dashed line) and are well separated from the softest 
magnons (states on the dotted line), (b) Spectrum of the Heisenberg model on the checkerboard 
lattice, which displays VBC LRO: for A'^ — > oo, the two lowest 5 = states (crosses) become 
degenerate and the gap to the others states remains finite, the finite spin-gap implies that SU (2) 
remains unbroken, (c) Spectrum of the MSE model (eq. k| with J2 = —2, J4 = 1) on the triangular 
lattice, a Type I Spin Liquid: four S = states (crosses) become degenerate as N ^ 00, the 
excitations above the degenerate ground-state remain gapped, (d) Spectrum of the Heisenberg 
model on the kagome lattice, a Type H Spin Liquid: there is a spin-gap and a gapless continuum 
of singlet states. 



the investigation of properties out of range of the ED approach-, 
approaches, however, is Hmited by the well-known sign problemlij''^ 



Monte-Carlo 
in frustrated 
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systems like those considered below. 

A similar analysis of the ED spectra has shown JJiat the Heisenberg model on the 
triangular lattice displays noncoUinear Neel LRO.ErH It also rejziealed occurences of 
"order by disorder" in the J1 — J2 model on the triangular latticeEland the J1 — J2—J3 
model on the honeycomb latticeJj In these situations the classical ground-state is 
a degenerate manifold of different kinds of 4-sublattice Neel order. Evidence of the 
selection of a collinear LRO by quantum fluctuations appears when increasing N, 
with the separation of the QDJS specific of collinear LRO from the states associated 
with 4-sublattice order. _ 

Calculations of spin-spin correlations, order parameters, spin susceptibilities,a 
stiffnessealj yield results consistent with the spectrum analysis but their finite scal- 
ing behavior is less favourable, (i.e. ~ 1/VN in Neel systems) and extrapolation 
to — > 00 could be rather inaccurate. Agreement between the scaling behaviors 
of ED and spin wave results can nevertheless help to circomvent this limitation. 



3 VBC LRO 

This class includes systems with LRO of singlet objects such as singlet bonds 
(dimers) or more extented ones (4-site plaquettes...). This kind of LRO usually 
breaks the space symmetry, so the spectrum must display a corresponding degen- 
eracy in the thermodynamic limit, but it does not break SU{2) which manifest by 
a finite spin-gap. Such features are visible in Fig nh_which displays the spectrum 
of the Heisenberg model on the checkerboard lattice:E3 two S — states appear in 
the bottom of the spectrum well separated from all others. A finite size analysis 
confirm the degeneracy of the two states and that the singlet and the spin gaps 
remain finite for N —^ 00. In the present case, symmetry properties of the two 
states indicate unambigously LRO of plaquettes on the void squares of the checker- 
board lattice, which is confirmed by calculations of correlation functions. Since 
the excitations are separated from the degenerate ground-state by a gap A, the 
thermal behavior of the spin susceptibility x and heat capacity C„ is exponentially 
activated ~ exp(— A/fcT) at low T. Due to the gap A, the ground-state energy in 
the different spin sector S evolves as E{S) Eq + AS + ■ ■ ■ at low S. This leads, 
in a continous Landau description to /(m) w /(O) + Am/2 — mh/2 + ■ ■ ■ in the 
presence of a magnetic field h: a finite field he — A is needed to induce a transition 
to states with a non zero magnetization, the magnetization curve m{h) starts with 
a plateau at zero m. 

Much attention has been given to the Ji — J2 model on the square lattice which 
for J2/J1 ^ 0.5 hasj-heeii suspected to display dimer or plaquette LRO, but the issue 
is still debated (seetaEJ and references therein). Besides the checkerboardJattice,lla 
we have found a VBC phase in the J1 — J2 model on the honeycomb latticcH In these 
three systems the VBC phase is located in the phase diagram near a collinear Neel 
phase. VBC LRO appears as a generic situation to be expected after destabilization 
by quantum fluctuations of a collinear Neel phase with increased frustration. This 
is in agreement with various analytical calculations, based on large-A^ SU{N) or 
Sp{N) expansiDna^aL non linear sigma-model theory taking into account the Berry 
phases effect sJi3'ESEZl But ED calculations, limited to small systems, does not allow 
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to investigate the phase transition between Neel and VBC LRO (see RefsllZl'tj for 
recently proposed scenarios). Experimental evidence of VBC LRXi has been found 
in the 2d insulating compounds: SrCu2{B03)2cli and CaViOg.c3 These systems 
are however somewhat special. The 2d lattices of magnetic ions in these compounds 
have lower symmetries than the above-mentioned lattices. The hamiltonian already 
favors the formation of singlet objects. There is no further symmetry breaking. 



4 Type I Spin Liquid 

The spectrum of type I Spin Liquid (see Fig [^c) seems at first sight very similar to 
the spectrum of a VBC. It has a a finite gap above a degenerate ground-state. So x 
and Cv are thermally activated and SU (2) is not broken. But we are lead to exclude 
VBC LRO because the few 5 = levels collapsing together to the ground-state do 
not have symmetries that could be explained in the framework a VBC picture 
and all correlation functions are found short-ranged. Indeed, in the two systems 
where we think that±ype I Spin Liquid behavior occurs: for the MSE model on th£ 
triangular latticfj'CJ and in the Ji — J2 — J3 model on the honeycomb lattice,El 
VBC LRO would require a larger degeneracy of the ground-state than the one 
observed. For instance, the simplest VBC LRO on the triangular lattice is 12-fold 
degenerate whereas the observed degeneracy is 4-fold. This 4-fold degeneracy isjJUt 
the signature of LRO in a local order parameter but has a topological origin.EffEj 
Several aspects of the physics of type I Spin Liquids, still deserve further inves- 
tigations. This fully gapped phase without VBC LRO is likely to be described with 
short ranged resonating valence bonds (RVB). A possible connection with the BVB 
phase of quantum dimer model, such as the one found on the triangular latticecJ is 
worth to be clarified. In XyvG I Spin Liquids, the excitations are expected to consist 
of deconfined spinous .Illrll3 We are presently searching for numerical evidence of 
deconfined spinous. It also remains to firmly outline the conditions of appearance 
of type I Spin Liquids. Fully gappejd-Spin Liquids are predicted to occur after 
destabilization of non-coUinear LRO.I13'll3 The two examples of type I Spin Liq- 
uid, found up to now, do not disagree with this. But, as discussed below, type 
II Spio-Liquid behavior with gapless excitations (not described in the approach of 
ReffO'Ea), may appear after destabilization of non-coUinear LRO. A further clue 
might be of importance: in the two type I Spin Liquids short range correlations 
are slightly ferromagnetic. In fact, they appears in the vicinity of a ferromagnetic 
phase. 

Experimentally, this phase might be observed in the low density 2d layers of 
solid ■^iJe where there is evidence that the MSE model on the triangulBX-Lattice 
with coupling constants that might be in the good range, isLrealized (seecj'OL^Ej 
and references therein) and possibly in the Wigner crystal.tZI It might also occur in 
the BaCo2 (^504)2 compound,Ej likely to be described by the Ji — J2 — J3 model on 
the honeycomb lattice with parameters close to those for which we found indications 
of this behavior. 
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Figure 2. Heisenberg model on the kagome lattice. Logarithm of the number of states in the 
spin-gap (squares fitted by the full line) vs size A'^ of the samples. The observed behavior differs 
from the scaling law exp(bN"^ ("+2) ^ deduced from a single mode description of the continuum 
with a dispersion law e(fc) = fc" for n = 1 (short-dashed line) and n = 2 (long-dashed line). The 
linear behavior observed here implies density of states increasing as (...)^ and a T=0 residual 
entropy. 



5 Type II Spin Liquid 

This_ahpiae Jifas first observed in the Heisenberg model on the kagome \at- 
tice.E3'OL3'Ej Some possible physical realizatians are SrCrGaO or various jarosites 
which have kagome planes of magnetic ions.EirE^ Fig|^d shows the spectrum of the 
A'^ = 36 sample of the kagome antiferromagnet. Two features, confirmed by finite 
size scaling, are visible. First, a finite, although small, spin-gap A. So SU{2) spin 
symmetry remains unbroken and x is exponentially activated. Second, a gapless 
continuum of singlet levels filling the spin-gap. Futhermore, the number of singlet 
levels in the spin-gap increases exponentially with the size N of the sample as shown 
m Fig I All the correlations (spin-spin, dimer-dimer, chiral-chiral, ...) are very 
small and short range. l3'0 

The understanding of the low lying continuum of jSing]ets_is still incomplete 
but somej)rogresses have been done in the last yearsOc^OcSl Recent ED cal- 
culationsEII on a kagome N = 24 sample depletted of two sites suggest that the 
elementary excitations likely consist of deconfined spinons: as shown in Fig |^, the 
spin-gap does not depend neither on the presence of the two nonmagnetic holes nor 
on their separation which suggest that spin- 1/2 excitations are essentially free (un- 
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Figure 3. Spin-gap of kagome N = 2i sample depletted of two sites vs the euclidian distance d 
between the two empty sites. The d = value is the spin-gap of the undepleted sample. 
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could be related to a family of short range RVB 
On the other hand, they cannot be described as 
Goldstone modes associated to a quasi LRO of dimers. Assuming non-interacting 
modes with a dispersion law e(fc) = fc", it is easy to obtain the internal energy 
and entropy versus temperature and thus the number of states below a given en- 
ergy. The logarithm of this number of states scales as . Fixing the unknown 
constants from the measured value of the number of singlets in the spin-gap for 
N = 12, we see in Fig ^ that the mode description is an order of magnitude off at 
N ~ 36. The exponential increase in the number of eigenlevels in the finite energy 
interval below the spin-gap suggests a finite entropy per spin at T=0. The same 
exponential increase in the number of eigenlevels in each spin sector can be at the 
origin of the anonialous^assy behavior "without chemical disorder" seen in sev- 
eral kagoni£^compoundsJ13 on which dynamical mean-field theory approaches might 
shed light .c2l A low lying continuum of singlets is consistent with the observationcj 
that the low tenincrature specific heat of SrCrGaO is essentially insensitive to large 
magnetic fields..L3 This picture is also compatible with the muons expeciments of 
Uemura et ac3 and the elastic spin diffusion measurements of Lee et a/.E3 |_. 

Type II Spin-Liquid behavictunight also occur in 3d, in pyrochlore compoundaS 
and gadolinium gallium garnet .c3_L|i|kewise, 2d systems of coupled Id chains might 
exhibit Spin-Liquid like behavior .e3e2I Lately, neutron scattering evMence for de- 
confined spinous on an anisotropic triangular lattice has been found.E3 

For a long time, it was widely believed that type II Spin-Liquid behavior was 
associated to special lattices like the kagome or pyrochlore lattices, which at the 
classical level display an infinite local ground-state degeneracy. We have now several 
counterexamples. The checkerboard lattice which is a 2d analog of the pyrochlore 
lattice and exhibit a large similar degeneracy at the classical level, was found to dis- 
play VBC LRO £3 A kagome-like spectra is observed on the triangular lattice when 
the 3-sublattice Neel order is destabilized by a frustrating 4-spin exchanged and we 
have no evidence of an infinite local degeneracy at this point. We now conjecture 
that the "kagome-like" phase is one of the possible scenario after destabilization 
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a non-coUinear LRO. At this point we have not made any difference between the 
criteria of appearance of the two Spin-Liquids. Attempting, to go further, one could 
speculate that type I Spin-Liquid behavior might originate from resonating entities 
including near neighbor triplet pairs assembled in singlet entities of at least 4 spins, 
whereas type II Spin-Liquid involve essentially near neighbor singlet pairing. A 
last hint comes from the study of the spin-S* = 1 kagome magnetjlj for sizes up to 
TV = 18, which point to a rather large spin-gap void of singlets. It is in fact a com- 
mon general idea that, due to topological reasons, half-integer and and integer spin 
systems could differ. The difference between spin-1/2 and spin-1 kagome spectra 
might possibly originate from this topological property. 
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